home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Libris Britannia 4
/
science library(b).zip
/
science library(b)
/
MEDICAL
/
H121A.ZIP
/
FILES6.EXE
/
lha
/
CENTILE.PGM
< prev
next >
Wrap
Text File
|
1991-07-15
|
3KB
|
53 lines
************************************************************************
*CENTILE.PGM Version 1.0, 20/FEB/1990 *
*An Epi Info program to present anthropometric centile distributions *
* by Kevin Sullivan, Division of Nutrition, CDC. *
* The Epi Info file analyzed must contain the anthropometry fields *
* HAZ, HAP, HAM, WHZ, WHP, WHM, WAZ, WAP, WAM, and FLAG. *
************************************************************************
echo
echo
echo You must read a file in order to run this program. Type
echo READ and the name of a file. There is an example data
echo file called NUTRI that can be used.
echo
?EPI>?
************************************************************************
*Declare missing variables based on Record Flag (same as AMISS1.PGM) *
************************************************************************
Define NOHA #
Define NOWH #
Define NOWA #
If (flag=1) or (flag=3) or (flag=5) or (flag=7) then NOHA = 1 else NOHA = 0
If NOHA = 1 then HAZ=.
If NOHA = 1 then HAP=.
If NOHA = 1 then HAM=.
If (flag=2) or (flag=3) or (flag=6) or (flag=7) then NOWH = 1 else NOWH = 0
If NOWH then WHZ=.
If NOWH then WHP=.
if NOWH then WHM=.
If (flag>3) and (flag<8) then NOWA = 1 else NOWA = 0
If NOWA then WAZ=.
If NOWA then WAP=.
If NOWA then WAM=.
************************************************************************
*Define field names for centile categories *
************************************************************************
Define HAPGRP _
Define WHPGRP _
Define WAPGRP _
************************************************************************
*Create centile categories using recode command *
************************************************************************
recode hap to hapgrp 0-9.99=a 10-19.99=b 20-29.99=c 30-39.99=d 40-49.99=e 50-59.99=f 60-69.99=g 70-79.99=h 80-89.99=i 90-99.8=j else= .
recode whp to whpgrp 0-9.99=a 10-19.99=b 20-29.99=c 30-39.99=d 40-49.99=e 50-59.99=f 60-69.99=g 70-79.99=h 80-89.99=i 90-99.8=j else= .
recode wap to wapgrp 0-9.99=a 10-19.99=b 20-29.99=c 30-39.99=d 40-49.99=e 50-59.99=f 60-69.99=g 70-79.99=h 80-89.99=i 90-99.8=j else= .
************************************************************************
*Ignore missing values for report *
************************************************************************
set ignore=on
************************************************************************
*Use report for centile output *
************************************************************************
report centile